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We report a study of the suppressed decays B~ [K'^ 7v~]dK~ and B~ — > [K'^'k~]dt'~ , where 
[K^n~]D indicates that the K^n~ pair originates from a neutral D meson. These decay modes 
are sensitive to the Unitarity Triangle angle (j)^. We use a data sample containing 275 million BB 
pairs recorded at the T(45') resonance with the Belle detector at the KEKB asymmetric e'^e~ 
storage ring. The signal for B~ — > [K^tv~]dK~ is not statistically significant, and we set a 
limit rs < 0.27 at 90% confidence level, where rs is the magnitude of the ratio of amplitudes 
\A{B^ D'^K^)/A{B^ D^K~)\. We observe a signal with 6.4a" statistical significance in the 
related mode, B^ [K^n^]]jT^^ ■ 

PACS numbers: 11.30.Er, 12.15.Hh, 13.25.Hw, 14.40.Nd 



Precise measurements of the elements of the Cabibbo- 
Kobayashi-Maskawa matrix Q constrain the standard 
model and may reveal new physics. However, the extrac- 
tion of the Unitarity Triangle angle (j)^ is a challenging 
measurement even with modern high luminosity B facto- 
ries. Several methods for measuring 03 use the interfer- 
ence between i?^ — > D'^K~ and D^K~ , which oc- 
curs when and _D" decay to common final states 0,0- 
As noted by Atwood, Dunietz and Soni (ADS) CP 
violation effects are enhanced if the final state is chosen 
so that the interfering amplitudes have comparable mag- 
nitudes; the archetype uses B^ [K^tt~]dK~ , where 
[if+TT^]/) indicates that the K'^tt^ pair originates from 
a neutral D meson. In this case, the color-allowed B de- 
cay followed by the doubly Cabibbo-suppressed D decay 
interferes with the color-suppressed B decay followed by 
the Cabibbo-allowed D decay (Fig. ^| . Previous studies 
of this decay mode have not found any signals |^. For 
the suppressed decay B~ [X+Tr^J^iTr", both topology 
and phenomenology are similar to B^ — > [K~^Tr~]]jK^ . 

In this analysis, the favored decays B^ 
[K^T:^]£)h^ , where h = n or K, are used as control sam- 
ples to reduce systematic uncertainties. The same selec- 
tion criteria for the suppressed decay modes are applied 
to the control samples whenever possible. Throughout 
this report, charge conjugate states are implied except 
where explicitly mentioned and we denote the analyzed 




FIG. 1: B' [K+tv']dK- and ^ [K+tt-]d-!t- decays. 



decay modes as follows: 

Suppressed decay B^ [if +7r^]£i/i^ B^ Dsuph~ 
Favored decay B^ [if ^tt+J^i/i^ B^ — > Dfavh~ 

The results are based on a data sample containing 275 
million BB pairs, collected with the Belle detector at the 
KEKB asymmetric energy e+e~ collider operating at 
the T(4S') resonance. The Belle detector is a large-solid- 
angle magnetic spectrometer that consists of a silicon 
vertex detector (SVD), a 50-layer central drift chamber 
(CDC), an array of aerogel threshold Cerenkov counters 
(ACC), a barrel- like arrangement of time-of- flight scintil- 
lation counters (TOF), and an electromagnetic calorime- 
ter (ECL) comprised of CsI(Tl) crystals located inside a 
superconducting solenoid coil that provides a 1.5 T mag- 
netic field. An iron flux-return located outside of the 
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coil is instrumented to detect mesons and to identify 
muons (KLM). The detector is described in detail else- 
where J\- Two different inner detector configurations 
were used. For the first sample of 152 million BB pairs, 
a 2.0 cm radius beampipe and a 3-layer silicon vertex de- 
tector were used; for the latter 123 million BB pairs, a 
1.5 cm radius beampipe, a 4-layer silicon detector and a 
small-cell inner drift chamber were used (S| . 

Neutral D mesons are reconstructed by combining 
two oppositely charged tracks. For each track, informa- 
tion from ACC, TOF and specific ionization measure- 
ments from the CDC are used to determine a K/-k likeli- 
hood ratio P{K/n) = Ck/{Ck + Ct^), where £k and 
are kaon and pion likelihoods. We use the parti- 
cle identification requirements P{K/tt) > 0.8 for kaons 
and P{K/tt) < 0.2 for pions. These requirements se- 
lect kaons (pions) with momentum dependent efficien- 
cies of 80 — 95% (90 — 95%) and pion (kaon) misidenti- 
fication probabilities of 5 - 20% (15 - 20%). D candi- 
dates are required to have an invariant mass within ±2.5(7 
of the nominal D mass: 1.850 GeV/c^ < M{Kii) < 
1.879 GeV/c^. To improve the momentum determina- 
tion, tracks from the D candidate are refitted accord- 
ing to the nominal D mass hypothesis and the recon- 
structed vertex position (a mass-and-vertex-constrained 
fit). B mesons are reconstructed by combining D can- 
didates with primary charged hadron candidates. The 
signal is identified by two kinematic variables, the energy 
difference AE = Ed+E)^- — £^boam and the beam-energy- 
constrained mass Mbc — V ^boam ^ (PD + Ph- )^, whcrc 
Ed is the energy of the D candidate, E/^- is the energy of 
the h~ , and Sbeam is the beam energy, all evaluated in the 
center-of-mass (cm) frame; pn andp/j- are the momenta 
of the D and h~ in the cm frame. Event candidates are 
accepted if they have 5.2 GeV/c^ < Mbc < 5.3 GeV/c^ 
and I A£^| < 0.2 GeV. If there is more than one candidate 
in an event, we select the best candidate on the basis of a 
determined from the difference between the measured 
and nominal values of M{Ktt) and Afbc. 

To suppress the large background from the two-jet-like 
e+e" ^ q<j {<} = u,d,s,c) continuum processes, vari- 
ables that characterize the event topology are used. We 
construct a Fisher discriminant of Fox- Wolfram mo- 
ments called the Super- Fox- Wolfram (SEW) where 
the Fisher coefficients are optimized by maximizing the 
separation between BB events and continuum events. 
Furthermore, cos^b, the cosine of the angle in the cm 
system between the B flight direction with respect to 
the beam axis is also used to distinguish BB events 
from continuum events. These two independent vari- 
ables, SEW and cos 6b, are combined to form a likeli- 
hood ratio Tl = Csig/ (Csig + £cont), where £sig and /!cont 
are likelihoods calculated from the SEW and cos 9b dis- 
tributions of signal and continuum background events, 
respectively. We optimize the TZ requirement by maxi- 
mizing S'/V S + N, where S and N denote the expected 



numbers of signal and background events in the signal 
region. For B~ DsupK~ {B~ Dsup""") we require 
n > 0.85 (7^ > 0.75), which retains 44.8% (57.6%) of the 
signal events and removes 96.2% (93.2%) of the contin- 
uum background. 

For B^ DsupK , there can be a contribution from 
B- D^Ti-, D° K+K-, which has the same final 
state and can peak under the signal. In order to reject 
these events, we veto events that satisfy 1.843 GeV/c^ < 
M{KK) < 1.894 GeV/c^. The favored decay B- 
Dfa.vh~ can also cause a peaking background for the sup- 
pressed decay modes if both the kaon and pion from the 
ZJfav decay are misidentified. Therefore, we veto events 
for which the invariant mass of the Kit pair is inside 
the D mass window when the mass assignments are ex- 
changed. After applying this veto, the residual back- 
ground from Kit misidentification is found to be negligi- 
ble. The three-body charmless decay B^ K~^K~w~ 
can peak inside the signal regions of AE and Mbc for 
B~ DsupK^. This background is estimated from 
the AE distribution of events in a Z) mass sideband, 
defined as 1.637 GeV/c^ < M{KTr) < 1.836 GeV/c^ 
and 1.893 GeV/c^ < M{Ktt) < 2.093 GeV/c^. The 
estimated peaking background inside the AE signal re- 
gion is 1.7 ± 0.9 events, which we subtract from the ob- 
served B^ DsupK~ yield. As a check, we estimate 
the expected background level from the measured B^ 
K~^K~TT^ branching fraction Using this result, the 
expected background in our signal region is 2.1 ± 0.6 
events, assuming that the B^ K^tt^ yield is uni- 

formly distributed in phase space. For B^ -~^ Z3sup7r~, 
the peaking background estimated from the D mass side- 
band is consistent with zero. 

The signal yields are extracted using binned maximum 
likelihood fits to the AE distributions of events in the 
Mbc signal region, (5.27 GeV/c^ < Mbc < 5.29 GeV/c^). 
Backgrounds from decays such as B~ Dp~ and 
B~ D*Tr~ are distributed in the negative AE region 
and make a small contribution to the signal region. The 
shape of this BB background is modeled with a smoothed 
histogram obtained from generic Monte Carlo (MC) sam- 
ples. The continuum background populates the entire 
AE region. The shape of the continuum background 
is modeled with a linear function. The slope is deter- 
mined from the AE distribution of the Mbc sideband 
(5.20 GcV/c2 < Mbc < 5.26 GeV/c^). The AE fitting 
function is the sum of two Gaussians for the signal, the 
linear function for the continuum, and the smoothed his- 
togram for the BB background distribution. 

In the fit to the AE distribution of B^ ZJfavTr^, the 
free parameters are the position, width and area of the 
signal peak, and the normalizations of continuum and 
BB backgrounds. For the signal, the relative widths and 
areas of the two Gaussians are fixed from the signal MC. 
For the B~ Dfa.-vK~ fit, the position and width of the 
signal peak are fixed from the B~ — > DfavTr" fit results. 
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TABLE I: Summary of the results. For the — > DsupK^ signal yield, the peaking background contribution has been 
subtracted. The first two errors on the measured production branching fractions are statistical and systematic, respectively, 
and the third is due to the uncertainty in the B~ Dfav^" product branching fraction used for normalization. 



Mode 


Product branching 
fraction from [12] 


Efficiency 
(%) 


Signal Yield 


Statistical 
significance 


Measured product 
branching fraction 


Upper limit 
(90%C.L.) 


B- 

B~ -+ 
B- 




(6.9 ±0.7) X 10"'^ 
(1.4 ±0.2) X 10"^ 
(1.9 ±0.1) X 10"" 


12.9 ±0.2 
20.1 ±0.2 
12.9 ±0.2 
20.3 ±0.2 


28.5 _7 4 
376.0 +211 
8181.9 


2.3cr 
6.4a 


(3.2i^-^±0.2±0.5) X 10"' 
(6.6ti;?±0.4±0.3) X 10"^ 


6.3 X 10"' 



To fit the feed-across from B~ — > ZJfavTr" , we use a Gaus- 
sian shape where the left and right sides of the peak have 
different widths since the shift caused by wrong mass as- 
signment makes the shape asymmetric. The shape pa- 
rameters of this function are fixed to values determined 
by the fit to the B~ — > DfavTr" distribution using a kaon 
mass hypothesis for the prompt pion. The areas of signal 
and feed-across from B" Dtt^ , and the normaliza- 
tions of continuum and BB backgrounds are floated in 



the fit. For 



DsupK- 



and B -DgupTr , the sig- 



nal and BB background shapes are modeled using the fit 
results for the B~ DfavK~ and B~ — > Dfav7r~ modes, 
respectively. The area of the feed-across from _Dsup7r~ 
is estimated as the measured yield of _B" _Dsup7r" 
multiplied by the tt to K misidentification probability. 
However, the areas of the signal and the normaliza- 
tions of continuum and BB backgrounds are floated. 
The fit results are shown in Fig. [3 The numbers of 
events for B~ — > _Dsup^~ a-nd Dfi^^h" , and the sta- 
tistical significances of the B" D^^ph^ signals are 
given in Table ^ The statistical significance is defined 
as y'— 2 ln(£o / 'Cmax) , where £max is the maximum like- 
lihood in the AE fit and Co is the likelihood when the 
signal yield is constrained to be zero. The uncertainty 
in the peaking background contribution is taken into ac- 
count in the statistical significance calculation. 
The ratio of branching fractions is defined as 



Roh 



B{B- 



Dsuph 



N 



- / f -Dsup/i- 



where No^^^h {Not^^h) and eo.^^h- {^Dt^^h-) are the 
number of signal events and the reconstruction efficiency 
for the decay B" Dg^ph^ [B^ —f Dfavh^), and are 
given in Tabled 

The ratios Rjjh are calculated to be 



Rdk = (2.3^J:^(stat) ±0.1(syst)) X 10"^ 
Rdtt = (3.5iJ:JJ(stat)±0.2(syst)) X 10"^ 

Since the signal for _B" ^ DgnpK^ is not significant, we 
set an upper limit at the 90% confidence level (C.L.) of 
Rdk < 4.4 x 10"^, where we take the likelihood function 
as a single Gaussian with width given by the quadratic 
sum of statistical and systematic errors, and the area is 




-0.10 0.00 0.10 0.20 -0.20 -0.10 0,00 0.10 0.20 

AE (GeV) AE (GeV) 

FIG. 2; AE fit results for (a) B" ^ D^upK- , (b) B~ ^ 
fsupTr", (c) B~ — > D{s.vK~ , and (d) B~ Dfav^r". Charge 
conjugate modes are included in these plots. 



normalized in the physical region of positive branching 
fraction. 

Most of the systematic uncertainties from the de- 
tection efficiencies and the particle identification can- 
cel when taking the ratios, since the kinematics of the 
B" Dsuph~ and B~ — > Z?fav^~ processes are similar. 
The systematic errors are due to uncertainties in the yield 
extraction (4.7% — 5.4%) and the efficiency difference be- 
tween B- ^ Dsnph- and B" ^ Aav/i~ (1.3% - 1.7%). 
The uncertainties in the signal shapes and the qq back- 
ground shapes are determined by varying the shape of 
the fitting function by ±lcr. The uncertainties in the 
BB background shapes are determined by fitting the AE 
distribution in the region -0.07 GeV < AE < 0.20 GeV 
ignoring the BB background contributions. The uncer- 
tainties in the efficiency differences are determined us- 
ing signal MC. The total systematic error is the sum in 
quadrature of the above uncertainties. 

The product branching fractions for B" Dsnph~ are 
determined as 



BiB- ^ Dsnph-) - BiB- ^ Dt,,h-) x Rdh, 
and are given in Tabled A third uncertainty arises due 
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(d) 




DsupK^ is sub- 



-0.20 -0.10 0.00 0.10 0.20 -0.20 -0.10 0.00 0.10 0.20 

AE (GeV) AE (GeV) 

FIG. 3; AE fit results for (a) B' D^upK' , (b) B+ 
D,upK+, (c) B- AupTT", and (d) B+ -» Dsup7r+. 



TABLE II; Signal yields and partial rate asymmetries. 



Mode 


N{B-) 


N(B+) 






B D.upK 

B DsupTT 


8.2 

18.81^:^ 




0.88 
0.30 


-0.62 ± 0.06 
t°-J^±0.06 



to the error in the branching fraction of B~ —^ Df^^^h^ ^ 
which is taken from 12] . The uncertainties are statistics- 
dominated. For the B~ DsupK^ branching frac- 
tion, we set an upper limit at the 90% C.L. of B{B~ 



DsupK~) < 6.3 X 10~^. For B~ DsupTT~ , our mea- 
sured branching fraction is consistent with expectation 
neglecting the contribution from B^ D^Tr~ . 
The ratio Rdk is related to (f>3 by 



Rdk 

where 



'^TBTd cos 03 COS (5, 



A{B- D°K-) 



A{B- DOR") 



A{D^ ^ K-n+) 



5 = Sb + Sd, 



0.060 ±0.003, 



and Sb and Sd are the strong phase differences between 
the two B and D decay amplitudes, respectively. Using 
the above result, we obtain a limit on r^. The least 
restrictive limit is obtained allowing ±la variation on 
and assuming maximal interference (03 — Q°,d — 180° 
or 03 = 180°, S = 0°) and is found to be rs < 0.27. 

We search for partial rate asymmetries Adh in 
Dauph"^ decay, fitting the B^ and i?+ yields separately 
for each mode, where Adh is determined as 



A 



Dh 



B{B- ^ Dsuph-) - B{B+ -> Dsuph+) 
B{B- ^ Aup/i-) + B{B+ ^ 7^sup/i+) ■ 



The peaking background for B^ 
tracted assuming no CP asymmetry. The fit results are 
shown in Fig. Eland Table IHl We find 



Adk = 0.88 
Ad^ = 0.30 



^°:^^(stat)±0.06(syst), 
:[J:25(stat)±0.06(syst), 



where the systematic uncertainties arise from possible 
biases in the analysis algorithms (estimated from the 
B^ DfavTr^ control sample to be 2.5%); uncertainties 
in the extraction of the B~ and B^ yields (estimated by 
varying fitting parameters by ±lcr to be 4.9%); asym- 
metry in the particle identification efficiency of prompt 
kaons (estimated in [l^l to be 0.6%). We assume no CP 
asymmetry in the peaking background, and do not assign 
any systematic uncertainty from this source [l^ . 

In summary, using 275 million BB pairs collected with 
the Belle detector, we report studies of the suppressed 
decay B^ Dsuph+ {h = K,Tr). We observe B^ 
DsupT^~ for the first time, with a significance of 6Aa. 
The size of the signal is consistent with expectation based 
on measured branching fractions The significance 

for B^ DsupK~ is 2.3a and we set an upper limit 
on the ratio of B decay amplitudes tb < 0.27 at 90% 
confidence level. This result is consistent with previous 
searches 1^, and with the measurement of in the decay 



B~ 



DK-, D 



jl5| . which provides the 



most precise current determination of 03. 
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